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ABSTRACT 

Here we show an example of a young asteroid cluster located in a dynamically stable region, 
which was produced by partial disruption of a primitive body about 30 km in size. We estimate 
its age to be only 1.9 ± 0.3 Myr, thus its post-impact evolution should have been very limited. 
The large difference in size between the largest object and the other cluster members means 
that this was a cratering event. The parent body had a large orbital inclination, and was subject 
to collisions with typical impact speeds higher by a factor of 2 than in the most common 
situations encountered in the main belt. For the first time we have at disposal the observable 
outcome of a very recent event to study high-speed collisions involving primitive asteroids, 
providing very useful constraints to numerical simulations of these events and to laboratory 
experiments. 
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1 INTRODUCTION 

The asteroid population, being st eadily subject to a process 
of collisional evolution jPavis et al . 1989; Bott ke et all l2005l : 



iMorbidelli et alj|2009l : IXsphaud l2009h . provides excellent possi- 
bilities to study physics of collisional events. Asteroid families, 
which are believed t o originate from cat astrophic disruption of 
single parent bodies feappala et al.' '200 2|), are, alm ost one cen- 
, tury since the pioneering work by .Hiravamal ( Il918h . still an at- 
' tractive and challenging subject. They provide a key to our un- 
derstanding of the collisional history of the main as teroid belt 
teottke et alj2005l ; ICellino, DeU'Oro. & Tedescoll2009h . outcomes 
of disrupti on events over a size range inaccessible to laboratory ex- 
periments (iMichel. Benz. & Richardson .2003 ; Durda et al. .2007i ; 
lAsphaudlioioir clues on the mineralogical structure of their par- 
ent bodies teellino et al.'2002^, the role of space weathering effects 
jNes vomv et al. 2005; Vemazza et al. 2009) and to many other sub- 
jects. 

So far, ejecta from a few tens of large-scale collisions has 
been discovered across the main asteroid belt (e.g. [Zappala et al.l 

1 19951 ; iMothe-Diniz. Roig. & Carvandl2005l ; iNesvomv et al.l2005h . 
In terms of their estimated ages, most families identified so far are 
fairly old and have had enough time to evolve significantly since 
the epo ch of their formation a s a consequence of (i) chaotic dif- 
fusion jNesvomv et al.l l2002al ; iNovakovic. Tsiganis. & Knezevid 

|2010b'), (ii) semi-major axis d rift due to Yarko vsky effect 
(Farinella & Vokiouhlic kvl Il999l; iBottke et al.l l200l[). (iii) sec- 
ondary collisions jMarzari. Farinella. & DavisI 1999 ; Bottke etsll 
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|2005|) . (iv) non-destructive collisions jDeU'Oro & Cellind [2007h 
and/ or (v) diffusion due to close encounters with massive aster- 
oids jCamiba et al.l2003l ; [Novakovic. Tsiganis. & Knezevicl2010d ; 
iDelisle & Laskaill2012h . 

In this respect, little altered recently bom families may 
provide more direct information about the physics of break- 
up events. Evidence of recent collisions in the asteroid belt 
have been reported in the last decade and our knowledge 
about young asteroid families has been increased significantly 
^Nesvomv et al."2002b; 'Nesvo mv. Vokrouhlickv. & Bottkdl2006al ; 
Nesvomv & Vokrouhlickv 2006). Most of these groups are formed 
by asteroids belonging to the S taxonomic class. There are, how- 
ever, several important differences among the S and C-type as- 
teroids. The objects belonging to former class are thought to have 
experienced some thermal evolution since the time of their for- 
mation, and it is, for example, known that space weatheri ng pro- 
cesse s are different for these two classes of objects (e.g. iGaffevI 
I2OIOI) . Also, numerical simulations show that the outcomes of col- 
lision al events are dep endent on intemal structure of the parent 
body jjutzi et al] 12009 ). Because of these reasons it is necessary 
to identify also young C-class families in dynamically stable re- 
gions, because a few such groups are already known, but none of 
these is well suited to extract reliable enough information. Two 
C-type families, namely Veritas and Theobalda, about 8.3 and 
6.7 Myr old respectively, are both located in dynamically unsta- 
ble region (Nesvomv et al. 2003; Novakovic 2010a). Thus, despite 
their young ages, these families evolved significantly since post- 
impact situation. Mos t of the asteroids belonging to Beagle family 
dNesvomv et al.l2008h . which is probably less than 10 Myr old, are 
located in dynamically relatively stable region. However, this group 
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is embedded in the large Themis family making distinction between 
the real members of the group and background objects very diffi- 
cult. Finally, the youngest known group that might be formed by 
C-type asteroids is Emilkowalski cluster , which is only 220 ± 30 
kyr old jNesvomv & Vokrouhlickvll2006h . However, it seems to be 
rather an X- than C-type group because albedos of its members are 
much higher than expected for C-type objects. For example, geo- 
metric albedo of aste roid (14627) Emilkowalski is 0.2013±0.0170 
dMasiero et alj201lh . 

Thus, it is of extreme importance to identify young fami- 
lies, that belong to the most primitive C class, that do not suf- 
fer from above mentioned problems. We have found the first ex- 
a mple of this kind to be the Lorr e clus ter, recently discovered 
bv lNovakovic. Cellino. & Knezevi Jj). According to existing 
color data its largest member, (5438) Lorre, is a primitive car- 
bonaceous C-class asteroid, which may contain organic materials. 
Moreover, the members of this cluster are located in dynamical sta- 
ble regi on and very tightly packe d in the space of proper orbital el- 
ements dKnezevic & Milanil2003i ), suggesting a likely young age. 
Therefore, its post-impact evolution should have been very lim- 
ited. This makes it a very promising candidate for different possible 
studies. Two crucial prerequisites for these studies are an accurate 
identification of its members, and a reliable estimation of its age. 
These are the questions we address here. 



2 LORRE CLUSTER 
2.1 Membership 

A dynamical criterion for family membership is based on distances 
among the objects in the space of proper orbital elements: semi- 
major axis (Op), eccentricity (Cp), and inclination (ip). Usually, for 
this purpose the hierarchi cal clustering method (HCM ) and 'stan- 
dard' metric (d) are used l lZapDalaet aI]|l990l 1 19941) . This metric 
is defined as 



d = nap + 2(5ep)2 + 2(5sm(ip))2 



(1) 



where nap is the heliocentric velocity of an asteroid on a circular 
orbit having the semi-major axis Sp. Sap = ttp^ — ap^, Sep = ep^ — 
ep2 , and 5 sin(zp) = sin(jpj ) — sin(ip2 ), where the indexes (1) and 
(2) denote the two bodies under consideration. The HCM connects 
all objects whose mutual distances (expressed in meters per second) 

are below a threshold value (dc). 

Following the method described in Knezevic & Milani' t200(f) 
we calculated synthetic proper elements for 148 asteroids located 
in a region somewhat wider than that occupied by the cluster. This 
region covers the following ranges in the osculating orbital ele- 
ments: 2.738 < a < 2.758 au, 0.13 < e < 0.39 and 23 < z < 
31°. The number of asteroids includes numbered, multi- and single- 
opposition object^ found in the recent version of catalogs of os- 
culating elements retrieved from the AstDys web pageQ Then, we 
applied the HCM to this set of proper elements, and we analyzed 
the number of dynamically linked objects identified at different mu- 
tual distances (Fig.[T](. In particular, this was done by changing dc 
from 10 to 200 m s^^ at discrete steps of 10 m s~^. At the lowest 
tested value of dc = 10 m s^^ the HCM links 14 asteroids with 



^ Although the orbits of single-opposition objects are less reliably known, 
we used them as well in order to find as many cluster members as possible. 
^ Asteroids Dynamic Site: |http://hamilton.dm.unipi.it/astdys27| 
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Figure 1. Number of asteroids associated with Lorre cluster as a function 
of cut-off distance dc . Unusually tightly packed members and sharp distinc- 
tion of the cluster from background population are two main characteristics. 



Table 1. Proper orbital elements of the asteroids belonging to the Lorre 
cluster In columns are given: semi-major axis (Op), eccentricity (ep), sine 
of inclination (siniip)), mean motion (n), frequency of the longitude of 
perihelion ig) and frequency of the longitude of node (s). 



Asteroid 


ap [auj 




sin{ip) 


u[°/yrj 


9 ["/yr] 


s ["/yr] 


5438 


2.74732 


0.26290 


0.47230 


79.0466 


9.4486 


-49.7809 


208099 


2.74694 


0.26314 


0.47241 


79.0630 


9.4207 


-49.7557 


2001 RF42 


2.74427 


0.26321 


0.47176 


79.1781 


9.5380 


-49.7308 


2001 XFi67 


2.74718 


0.26314 


0.47253 


79.0532 


9.3965 


-49.7490 


2003 BW5 


2.74796 


0.26294 


0.47198 


79.0210 


9.4916 


-49.8129 


2003 YY120 


2.74671 


0.26342 


0.47212 


79.0732 


9.4846 


-49.8060 


2005 YDig 


2.74788 


0.26313 


0.47246 


79.0204 


9.4144 


-49.7840 


2006 AL 16 


2.74636 


0.26342 


0.47211 


79.0883 


9.4880 


-49.7946 


2006 RMgs 


2.74263 


0.26276 


0.47201 


79.2495 


9.4847 


-49.6157 


2007 BJ62 


2.74626 


0.26338 


0.47204 


79.0927 


9.4968 


-49.7919 


2008 AD104 


2.74722 


0.26308 


0.47240 


79.0511 


9.4276 


-49.7716 


2OIOCG1713 


2.74536 


0.26292 


0.47195 


79.1313 


9.4935 


-49.7350 


2OIIFQ151 


2.74521 


0.26299 


0.47196 


79.1377 


9.4915 


-49.7297 


2010 AX32 


2.74668 


0.26382 


0.47227 


79.0744 


9.4695 


-49.8253 


2006 VZ122 


2.74783 


0.26342 


0.47239 


79.0216 


9.4505 


-49.8372 


2008 BBio 


2.74631 


0.26348 


0.47220 


79.0899 


9.4758 


-49.7979 


2008 DEg 


2.74490 


0.26318 


0.47193 


79.1514 


9.5111 


-49.7482 


2010 EW42 


2.74544 


0.26333 


0.47198 


79.1277 


9.4987 


-49.7634 


2010 EJgi 


2.74233 


0.26344 


0.47210 


79.2627 


9.4890 


-49.6575 



Lorre, while the number of members raises to 19 for 20 m s ~^ . The 
number of dynamically associated members remains a constant un- 
til 80 m s~^, when one body, asteroid 2006 AXe? is added. Later 
on, no additional body is linked to the cluster, even for the largest 
used value of dc = 200 m s~^. 

From these results we can draw three basic conclusions: (i) 
the cluster is extremely compact and very well separated from the 
background population; (ii) the nominal membership of the cluster 
is best characterized at dc = 20 m s^^; (iii) the asteroid connected 
with the group at 80 m s^^ is likely a close background object. 
Thus, the Lorre cluster has 19 currently known members (Table[Tll. 
These asteroids are very tightly packed with mutual distances sig- 
nificantly smaller than in the cases of typical families in the main 
asteroid belt. 
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Table 2. Lyapunov times of Loire cluster members derived using different 
dynamical models. 
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Figure 2. The Lorre cluster in the space of proper elements. The size of each 
symbol is proportional to the diameter of the body. The superimposed el- 
lipses represent equivelocity curves, computed according to the Gaussian 
equations I Nesvomy & Vokrouhlickv 2006). These ellipses are obtained 
assuming a velocity change Ati = 15 m s~^, argument of peiihelion 
Lo = 90°, and true anomaly / = 90°. The ellipses are shown to illustrate 
distribution of the fraginents in the case of isotropic ejection field; however, 
it is easy to see that the ejection velocity field of the Lorre cluster was highly 
asymmetric, a nice example of what one should expe ct to be the outcome 
of a cratering event l lVokrouhlickv & Nesvomvl201ll) . The locations of the 
relevant mean motion resonances are denoted with the vertical dashed lines, 
except in the case of 3 J— IS— 1 A resonance that is marked with the gray- 
shaded region. In the top panel basic information about parent body are also 
given. 



2.2 Dynamics 

The orbits of the asteroids belonging to the Lorre cluster are char- 
acterized by the moderate eccentricities (ep « 0.26) and high in- 
clinations {ip ~ 28°), but the region occupied by these asteroids is 
not under influence of any of the strong mean motion or secular res- 
onances. Thus, despite their orbital characteristics, these asteroids 
are mostly stable. Still, there are a few mean motion resonances 
(MMRs), present in the region, whose influence should not be ne- 
glected. The most powerful is a three-bod)0 MMR 3 J- IS - 1 A lo- 



Three-body mean motion resonances are co mensurabilities between 
the m ean motions of Jupiter, Saturn and asteroid l lNesvomv & Morbidellil 
Il998h . They are characterized by the relation mjXj + ms^s + ~ 0> 
where Aj, A5 and A denote mean motions of Jupiter, Saturn, and asteroid 
respectively, while mj,mg and m are integers. 



Asteroid 


4pla 


7pla 


7 pla + Ceres 


7 pla + CV 


7 pla + CPV 


5438 


107.6 


41.7 


35.0 


27.2 


26.0 


208099 


128.1 


754 


50.6 


32.5 


30.7 


2001fl.F42 


290.7 


57.7 


36.0 


29.0 


16.1 


2001XFie7 


82.1 


48.2 


314 


35.1 


30.6 


2003 SVK5 


41.2 


6.7 


21.8 


7.9 


20.5 


2003^1-120 


534.8 


203.3 


39.9 


30.0 


32.7 


20051-018 


32.5 


38.9 


29.3 


25.1 


20.9 


2OO6AL16 


304.0 


162.1 


41.5 


35.8 


314 


2006 HM98 


226.2 


22.1 


14.0 


19.2 


19.7 


2007SJ62 


289.1 


219.3 


39.7 


34.9 


35.1 


2008AJ3io4 


125.6 


41.6 


37.8 


28.6 


26.0 


2010CGi76 


7142.9 


75.1 


374 


36.9 


23.0 


2OIIFQ151 


505.5 


115.1 


34.4 


27.1 


24.7 


2OIOAX32 


1960.0 


198.4 


40.9 


35.2 


31.0 


2OO6V-Z122 


354 


31.8 


27.5 


25.1 


il.5 


2OO8BB10 


552.9 


202.0 


43.6 


35.6 


29.7 


200813 Fg 


4000.5 


83.8 


37.6 


35.3 


31.3 


2OIOFIV42 


2381.0 


76.0 


44.6 


337 


28.8 


2OIOFJ81 


746.6 


97.9 


42.3 


34.2 


34.1 



cated at 2.752 au. A somewhat less significant but still relevant are 
another two 3-body MMRs, namely 1J-I-4S-1A and 4J-I-3S-2A 
(see Fig.[2ll. Finally, 13J/5A 2-body MMR, among Jupiter and as- 
teroid, is present in the region as well. 

To better understand the strength of these resonances and their 
possible influence on the dynamical stability we have determined 
Lyapunov times (Tiyap) for all members of the L orre cluster. This 
was d one according to the method proposed by iMilani & Nobilil 
and within the framework of several different dynamical 

models. 

As for most of the purposes, in this part of the main aster- 
oid belt, dynamical model with four major planets (from Jupiter 
to Neptune) is accurate enough, we first used this model to esti- 
mate Lyapunov times. The obtained values of Tiyap are in most 
cases longer than 100 kyr. A few exceptions include objects lo- 
cated around ap — 2.7478 au, that are probably trapped inside the 
4J-H3S— 2A resonance. However, even Lyapunov times of these ob- 
jects are not shorter than about 30 kyr (Table|2}- 

When dynamical model with seven planets, from Venus to 
Neptune, is used, the estimated Lyapunov times are noticeably 
shorter (Table |2j, meaning that this model should be used for as- 
teroids located in the region of Lorre cluster The reasons for the 
important difference among the results obtained with 4- and 7- 
planets are relatively large orbital eccentricities and inclinations of 
these objects. Still, according to this result most of the Lorre cluster 
members are reasonably stable, with the only one possible excep- 
tion, asteroid 2 003 BW5. 

Recently, iLaskar et al.l ilOl ll) showed that close encounters 
with massive asteroids may induce chaos in their and in the mo- 
tion of other asteroids. To check whether or not this is the case for 
Lorre cluster members, we have also calculated Lyapunov times 
using dynamical models that include some of the most massive as- 
teroids, Ceres, Pallas and VestaQ 

Our result generally confirms that obtained by iLaskar et al.l 

* For this purpose, the masses of Ceres, Vesta and Pallas 
are set to 4.757, 1.300 and 1.010 X 10~^°Mfn respectively 
iKuzmanoski, Apostolovska. & Novakovid[2010l ; lBaer, Cheslev, & MatsonI 
l201lh . These masses are results of the latest calculations performed by 
means of the improved methodology. A preliminary estimation of Vesta's 
mass provided by Dawn mission (http://dawn.jpl.nasa.g0v/mission/l per- 
fectly match the results from these two papers, for this object. Due to these 
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ij). Lyapunov times become, on average, shorter when the mas- 
sive asteroids are included in the dynamical model. There are, how- 
ever, a few asteroids whose motion seem to be more stable in this 
case, and their values of Tiyap are longer, than those obtained in the 
model with 7-planets only. An illustrative example is the only pos- 
sibly unstable object among the currently known members of the 
cluster, the asteroid 2003 BW5. Its estimated Tiyap is only 7 kyr in 
the dynamical model with 7 planets, but rises to 22 kyr when Ceres 
is added to the dynamical modellf] Thus, although influence of the 
massive asteroids on the motion of asteroids belonging to the Lorre 
cluster is undoubtedly confiiTned, its resulting effect may vary from 
case to case. 

The conclusion that we can draw from derived values of Lya- 
punov times is that orbits of the Lorre cluster members are neither 
perfectly stable nor strongly chaotic. 

In terms of a possible post-impact dynamical evolution of the 
cluster even week chaos may be important. Hence, to explore this 
possibility and to assess a jet like shape of the cluster, we checked 
stability of the proper eccentricity and inclination of asteroids be- 
longing to the Lorre cluster. Using the numerical integrations of 
cluster members performed in the dynamical model that includes 
seven planets (from Venus to Neptune) and three most massive 
asteroids (Ceres, Pallas and Vesta), we estimated average evolu- 
tion rates of eccentricity and sine of inclinations to be 1 x 10^* 
and 5 x 10""" per one million years respectively. These are slow 
changes that do not seem to be able to significantly change over- 
all structure of the cluster. Actually, as we found the Lorre cluster 
to be only about 1.9 Myr old (see Section [Z3t , over its lifetime ex- 
pected changes of eccentricity and sine of inclination are only about 
2 X 10"* and 1 x 10"'' respectively. By comparing these values 
with the scales of y-axes in Fig. |2] we concluded that dynamical 
evolution is negligible. 

Looking at Fig.|2]it can be easily realized that distribution of 
the Lorre cluster members is highly asymmetric with respect to the 
largest member, asteroid (5438) Lorre. To understand the reasons 
for this, we extend our dynamical analysis to the region surround- 
ing the cluster. The dynamical instability starts to increase for val- 
ues of semi-major axis larger than 2.748 au. The inner border of the 
powerful 3J— IS — lA MMR is found at about 2.749 au. However, 
using numerical integrations of 100 massless test particles we have 
verified that this instability cannot explain the absence of cluster 
members in the 2.748 - 2.754 au range (see Fig.|2}. Although, over 
a time scale of 2 Myr, many particles interact with the 3J— IS— 1 A 
resonance, they still remain close enough to be recognized by the 
HCM. 

Available evidence suggests therefore that the observed asym- 
metry of the family is mostly a consequence of the original ejection 
velocity field of the fragments, rather than dynamical post-impact 
evolution. Thus, this cluster still keeps memoiy of the original ejec- 
tion velocity field, a useful input to study impact physics. 

2.3 Age 

The most accurate method known so far to estimate the age of 
a young asteroid family is to integrate the orbits of its mem- 
reasons we c hose to use these va lues, despite being slightly smaller than 
those used bv lLaskaret"aIU201lh . 

^ This is not a surprise because an estimation of Lyapunov times, even 
for moderately chaotic orbits, is probabilistic, thus not highl y reliable and 
should be interpreted with a carel lKnezevic & Ninkovidl2005l) . 



hers back wards in time and to i d entify the epoch of their con- 
vergence jNesvom^ etal.l l2002bl |2003|) . However, this method 
can be applied only to the objects on stable orbits. As we 
showed in Section 12.21 that the orbits of the Lorre cluster mem- 
bers are not perfectly stable an application of the backward in- 
tegration method (BIM) is not so straightforward. To overcome 
this problem we turn to a statistical approach based on the 
BIM jNesvomv & Vokrouh lickv 2006: Vo krouhlickv & Nesvomvl 
I20T 11). Instead of orbits of nominal members we used a number of 
cloned, statistically equivalent, orbits. In this way we were able to 
characterize the age of the Lorre cluster in a statistical sense. 

More in particular, we took into account the current orbital un- 
certainties of the nominal orbits and different possible evolutions of 
the orbital semi-major axes due to the Yarkovsky effect. For each 
nominal member of the cluster, except for asteroid (5438) Lorre, 
we produced a set of 10 orbital clones. These clones are drawn 
from 3a interval of their formal uncertainties listed in Table|3] as- 
suming Gaussian distribution. Then, for each of the orbit clones we 
generated 10 different 'yarko' clones uniformly distributed over the 
interval stretching fr om zero to the max imum expected drift due to 
the Yarkovsky force ( Bottke et al .120011) . The maximum drift in the 
proper semi-major due to the Yarkovsky force (da / dt)max for each 
object is obtai ned assuming thermal parame ters appropriate for C- 
type asteroids jBroz & Vokrouhlickvll2008l) . In this way a total of 
100 statistically equivalent clones were assigned to each member. 
Clones are not used for asteroid (5438) Lorre itself because on one 
hand its orbit is very well determined, while on the other hand it is 
large enough (see Table |4} that Yarkovsky effect on its orbit can be 
safely neglected. 

The orbits of all clones were numerically integrated backward 
in time for 10 Myr using the Orbit9 software. These integrations 
were performed within the framework of a dynamical model that 
includes seven planets, from Venus to Neptune, as perturbing bod- 
ies, and accounts also for the Yarkovsky effect^To account for the 
indirect effect of Mercury, its mass is added to the mass of the Sun 
and the barycentric correction is applied to the initial conditions. 

The age of the cluster was estimated by randomly select- 
ing one clone for each member and determining the age for that 
particular combination of clones as the minimum of the function 
(Vokrouhlickv & Nesvomy 2011): 

AV = nay^{sin{i)Any + 0.5{eAzuy (2) 

where na ~ 18 km s"' is the mean orbital speed of the asteroids 
in the cluster, and AQ and Azu are the dispersions of the longitude 
of node and the longitude of perihelion, respectively. 

The obtained results are shown in Fig.|3] The age of the Lorre 
cluster turns out to be 1.9 ± 0.3 Myr. The estimated error comes 
mainly from the assumed orbital uncertainties of single-opposition 
asteroids. Nevertheless, the result is robust and undoubtedly con- 
firms that the LoiTe cluster is very recent. 



^ For single-opposition objects we used the following values for all objects: 
Ua = 2.0 X IQ-^ au, o-e = 3.0 X IQ-^, <ji = 1°0 X IQ— an = 
2?0 X lO""*, (7u, = 3?5 X 10"* and ctm = 5°0 x IQ-^. 
^ For simplicity, the Yarkovsky effect is included in the model as a constant 
seculai' drift (inwai'ds or outwards) of the semi-major axis. This approxi- 
mation seems appropriate for our purpose to characterize the age of Lorre 
cluster in a statistical sense. 
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Table 3. The osculating orbital elements along with their formal uncertainties of Lorre cluster members at epoch 56000.0 MJD as found at AstDys. 
The horizontal line separates single-opposition from multi-opposition and numbered asteroids. 



Asteroid 


a [au] 


e 








M[°] 




a a 










M 


(5438) Lorre 


2.7457268384 


0.2763423275 


26.57394988 


29851646410 


238.57467176 


53.56378415 




0.0000000219 


0.0000001256 


0.00001231 


0.00001893 


0.00002667 


0.00002319 


(208099) 2OOOAO201 


2.7473300263 


0.3276837375 


24.27112923 
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Table 4. Different characteristics of the LoiTe cluster members. 
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Figure 3. The histogram of possible ages of Lorre cluster. It is constructed 
using 10^ different combinations of clones (see text). The age of cluster 
derived from these values is 1.9 it 0.3 Myr. 



2.4 Physical and spectral characteristics 

As for physical properties, the geometric albedo s {pv) have been 
determined for 6 m embers of the Lorre cluster dUsui et aTlbOl ll : 
iMasiero et ani2oT l|), with an average value of 0.053, compatible 
with C-class objects. 

Unfortunately, little is known about the spectral reflectance 
properties. To date, a spectral class has been determined only for 
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the largest asteroid (5438) Lorre, which is classified as a C-type 
( IBus & Binzell2002l) . 

It is interesting to note that for this same asteroid an estimate 
of the rotational period P is also available. According to ' Behreridl 
(201 1) P is about 25 hours. This unusually long period might 
be, at least part ly, the result of angular momentum transfer dur- 
ing the impact jPobrovolskis & Bumslll984l : ICellino et ^Il990l : 
iTakeda & Ohtsukill2009h . that may produce in some cases despin- 
ning of mid-sized objects. 
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2.4.1 Size of the parent body 

To further characterize the event which produced the Lorre clus- 
ter we estimated the size of the parent body. The simplest way to 
achieve this goal is to estimate the volume of the parent body by 
summing up the volumes of all known members, assuming a spher- 
ical shape for all of them. For this purpose we used the available 
diameters of the objects obtained by thermal radiometry observa- 
tions, using WISE data in all cases given in Table [T] but for Lorre 
itself, whose diameter is known from AKARI observations. For the 
objects lacking a size estimate (the majority of the objects in our 
sample), we derived it using the well known relation between di- 
ameter, absolute magnitude and albedo (see below). 

We adopted for each object the nominal value of its absolute 
magnitude H taken from the AstDys catalog, (these data are also 
listed in Table [T] One should be aware that the catalog values of 
H are known to be affected by large unce rtainties for objects in 
this magnitude range jMuinonen et al.l2010h . This also affects neg- 
atively the errors in the albedo determined by means of the thermal 
radiometry technique, and for this reason we tend to believe that 
the nominal values listed in Table[T]for WISE and AKARI-derived 
albedos may well be quite optimistic in some cases. For each object 
lacking an albedo measurement, we adopted the average value of 
0.053 for this family, which is based on the nominal values shown 
in Table[TJ. From H and the albedo we can derive the size from the 
relation log(L>) = 3.1236 - 0.2H - 0.5 log(pv) where D is the 
diameter. The obtained D values range between 1.1 and 4.1 km. 

By summing up all the resulting volumes of the family mem- 
bers, we find that the parent body was just a little larger than the 
largest fragmenjf], (5438) Lorre, which has an estimated diameter 
about 30 km. This conclusion does not change if we simply assume 
that the parent body could not be smaller than the sum of the sizes 
of the two largest family members. This is the criterion applied by 
lTangaetal.l ( ll999l) . and it is based on simple geometric considera- 
tions, which is more suitable to treat the cases of full parent body 
disruption. 

The escape velocity from a surface of a 30-km body is about 
13.5 m s^^ (assuming a density of 1.5 g cm^'^, typical of C- 
class asteroids). The second largest member of the cluster, asteroid 
(208099) 2000 AO201 . is about 6 km in diameter. The cluster turns 
out to be therefore the outcome of a cratering event, which was not 
sufficiently energetic to completely disrupt the parent body. This 
result supports our conclusion that the observed asymmetry of the 
cluster is likely a consequence of the original ejection velocity field. 



2.4.2 Ejection velocity field 

The structure of the families in the space of proper elements can 
be used to infer some information on the ejection velociti es of the 
fragments in family-forming events ('Zappala et al.ll2002h . As we 
already noted, the most important feature of the ejection velocity 
field (EVE) of Lorre seems to be a high asymmetry with respect to 
the location of the largest member. This, however, is not the only 
peculiar characteristics of the EVE. A jet like structure is visible in 
both, (ap,ep) and (ap,sin{ip)), planes (Fig.|2}. This is not unex- 
pected in the case of a cratering event. Jetting is expected to have a 
chance to occur when two objects collide at high speeds and at high 



° The difference among diameters of the parent body and largest fragment 
is smaller than the uncertainties of these two values. 



Table 5. Differences in velocities with respect to asteroid (5438) Lorre. 
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incidence angles (see e.g. iHousen & HolsappldlioTll) . Such struc- 
ture is obs erved in both numeric al simulations and laboratory ex- 
periments jYang & Ahrens|[T995h . but it has not been observed yet 
among real asteroid families, mainly due to the post-impact evolu- 
tion of the known groups. 

Although a detail study of the EVE is beyond the scope of this 
paper, we want to emphasize here that there is a clear trend in the 
velocity -size rela tionship. This trend is in agreement with previous 
studies jCellino et al. 1999) suggesting that smaller fragments are 
ejected on the average with slightly higher velocities. However, the 
number of known cluster members is still too small at the moment 
to analyze this trend in more detail. 

Finally, it should be noted that differences in velocities (Ta- 
ble O are much smaller than what is usually expected in the 
cases of dynamical families produced by disruption events. In fact, 
Lorre seems to be likely issued from a moderate-energy cratering 
event, and is the m ost compact group known so far amon g high- 
inclination families jNovakovic, Cellino. & KnezevidboTlh . 

2.4.3 Size-frequency distribution 

Some important information about the impact physics can be ob- 
tained by study ing the size-frequency distributions ( SEDs) of as- 
teroid families ( Tanga et al. I I 19991 : 1 Durda et al. I l2007h . It is gener- 
ally found that these distribution can be described by a power law, 
N(>D) oc D^" . Younger asteroid families generally have steeper 
SEDs which are generally thought to evolve with time toward shal- 
lower trends due to coUisional and dynamical erosion of the family. 
A correct way to fit these distributions, i.e. to estimate exponent a, 
is to adopt an approach based on maximum li kelihood method ap- 
plied to bi-trunca ted Pareto distributions (see ICellino et alj|l99ll : 
lTangaetalJi99^ . However, the number of family members is cur- 
rently too small to perform such a statistical analysis. 

Thus, we used an alternative approach based on the least- 
squares methocjflto estimate the exponent a. In this way, by fitting 
cumulative size distribution, for objects between 3.0 and 4.5 km 

^ This approach, despite being widely used, is not correct strictly speaking. 
In particular, this method may significantly underestimate the uncertainties 
of obtained values. However, as in any case the number of known members 
is too small to obtain a highly reliable result, we used this method because 
of simplicity of its implementation. 
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Diameter, D [km] 

Figure 4. Cumulative size distribution N(>D) of 19 asteroids belonging to 
the Loire cluster. 

in diameters, we found a to be 3.2. This v alue is smaller than ex- 
pected for typical young asteroid families dNesvomv et al.ll2006bl : 
IParker et al.ll2008b . Likely, this result is affected by the observa- 
tional incompleteness and a real a is somewhat larger. In any case, 
for a moment, we can only say qualitatively that the cumulative size 
distribution does not appear to be very steep. 

2.4.4 Collisional lifetime 

It is interesting to estimate what was expected collisional lifetime of 
the Lorre cluster parent body. This computation depends on many 
parameters, including mainly the inventory and size distribution of 
the possible impactors, the average impact velocity and, for what 
concerns the outcomes of the collisions, on the impact strength of 
the body, which in turn depends on its size and density. 

We computed the mean intrinsic collision probability and 
the mean impact velocity for the collisions between (5438) 
Lorre and other main-bel t asteroids using the approach of 
bell' Pro & PaoUcchil ^\99^ . The mean impact velocity results to 
be about 10 km s~^, due to the high-inclination orbit of (5438) 
Lorre. Under standard assumptions on the cumulative size dis- 
tribution of the population of possible projectiles, described by 
a power-law with an exponent of 2.5, a density value of 1.5 
g cm~ ^, and setting the impac t strength on the basis of the re- 
sults of lBenz & Asphaud ( Il999f) . the estimated collisional lifetime 
of (5438) Lorre is 6.6 Gyr, in ag reement with results of some inde- 
pendent studies teottke et all200 5). This relatively high value does 
not change much by steepening the size distribution of the projec- 
tiles (the lifetime becomes 5.3 Gyr if the power-law exponent is 
increased to the value of 3.0), nor by changing the value of the 
density. 

Asteroid Lorre is isolated, and there are no asteroids of sim- 
ilar size in its surroundings which might have been produced by 
the disruption of a hypothetical common parent body. We are led 
therefore to conclude that Lorre could be a pristine asteroid, which 
survived nearly intact since the time of its formation. This makes 
its analysis even more interesting. 



3 SUMMARY AND CONCLUSIONS 

Here we show the first example of a young asteroid cluster located 
in a dynamically stable region, which was produced by partial dis- 
ruption of a primitive body about 30 km in size. We estimate its 



age to be only 1.9 ± 0.3 Myr, thus its post-impact evolution is 
very limited. The large difference in size between the largest ob- 
ject and the other cluster members means that this was a crater- 
ing event. The parent body had a large orbital inclination, and was 
subject to collisions with typical impact speeds higher by a fac- 
tor of 2 than in the most common situations encountered in the 
main belt. For the first time we have at disposal the observable out- 
comes of a very recent event to study high-speed collisions involv- 
ing primitive asteroids, providing very usef ul constraints to numeri- 
cal simulations of these events |M ichel, Be nz. & Richardsorill2003l : 
lJutzi et al.l l2009':'Leinhardt & Stewart 2012) and to laboratory ex- 
periments IHousen & Holsapple 201 1). 

This is the best preserved young asteroid family produced by 
partial disruption of a primitive asteroid, of a kind which is sup- 
posed to have survived nearly unaltered since the epoch of forma- 
tion of the Solar System. Being young and well distinct from the 
background population, this cluster provides very useful informa- 
tion that can help to answer several long-debated questions in plan- 
etary science. Examples include a better understanding of impact 
physics, material strength and the role of space weathering. These 
process, highly dependent on the composition of the objects, are so 
far poorly constrained for primitive asteroids. 

Among the members of the Lorre cluster there are several as- 
teroid pairs, couples of objects with nearly identical orbital param- 
eters. These pairs may well consist of couples of fragments which 
were ejected with nearly identical ejection velocities. Another pos- 
sibility is that they might actually be the components of former 
binary systems originally produ ced by the collisio n, and later de- 
coupled by some mechanisms dPravec et al.ll201(il) . Production of 
binary systems i n collisional events has been suggested by numer- 
ical simulations ( lMicheletal1 l200l': Durd a et al.l 112004), but their 
expected abundance in asteroid families has not been firmly estab- 
lished yet. The young age of the Lorre cluster as well as its sharp 
separation from background objects may potentially help to better 
understand both populations, binaries and pairs. 

An interesting possibility for fu ture work comes from a re- 
cent result of iBenavidez et a l. ' 12012") how found that low-energy 
impacts into rubble-pile and monolithic targets produce different 
features in the resulting SFD, and, thus, this is a potentially di- 
agnostic tool to study the initial conditions just after the impact 
and the internal struct ure of th e paren t bodies of asteroid families. 
According to Benavidez et al.l ( |2012|) . cratering events, produced 
by small impactors, can potentially provide even more information 
about the internal structure of the parent body than catastrophic or 
super-catastrophic events produced by large impactors. Thus, the 
LoiTe cluster seems to be a very promising candidate. 

Next, the Lorre cluster may be very useful to improve our 
knowledge about space weathering processes acting on primitive 
bodies, a debated subject since results based on the Sloan D igi- 
tal Sky Survey broadband photometry dNesvomv et al.ll2005h are 
not consistent with the results of some laboratory experiments 
( lBnmettdl2009l) . 

Finally, the cluster may be a very interesting place to search 
for new main-belt comets (MBC s). F°l A recent findings by 
iNovakovic. Hsieh. & Cellindl20T^) supports an idea that this kind 
of objects may be pre ferentially found among the mem bers of 
young asteroid families dNesvomv et al.ll2008l:lHsiehll2009[) . In this 



Main belt comets are objects dynamically indistinguishable from main 
belt asteroids, but which exhibit com et-like activity due to the sublimation 
of volatile ice iHsieh & Jewitdl2006h . 
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respect, members of Lorre cluster are particularly interesting can- 
didates because their heliocentric distances are smaller than those 
of currently known MBCs. Thus, they may provide a clue about the 
inner edge of populations of MBCs. 
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